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In  the  phenomena  of  aircraft  flight,  the  body  is  subjected  to  forces 
and  stress  that  are  not  ordinarily  encountered.  Examples  are  emergency 
egress  from  aircraft,  high  G  turns,  and  the  random  vibrations  and  dis¬ 
turbances.  Anatomically  the  vertebral  colunn  must  absorb  much  of  these 
excesses.  In  light  of  these  biomechanics  problems  of  flight,  this  thesis 
will  study  the  influence  of  surgical  herniation  on  the  viscoelastic 
properties  of  the  intervertebral  joint.  The  presentation  shall  be  broken 
down  into  two  phases:  experimentation  and  analysis. 

The  experimentation  involved  several  distinct  phases.  First  five 
specimens  were  excised  from  the  L ^ - L ^  level  of  young  adult  rhesus 
monkeys  (Macaca  mulatta).  The  prepared  specimen  consisted  of  the  inter¬ 
vertebral  joint  and  the  two  adjacent  vertebral  centrums  without  either 
the  posterior  elements  or  associated  soft  tissues.  The  prepared  specimen 
was  subjected  to  a  constant  compressive  load  of  sixty-seven  newtons  for 
eight  hours  followed  by  a  relaxation  time  of  sixteen  hcu^s.  The 
specimens  were  then  subjected  to  various  stages  of  herniation  and  the 
creep  tests  were  repeated  after  each  herniation. 

The  analysis  was  also  carried  out  in  several  phases.  First  the 
experimental  data  was  graphed  using  a  specially  prepared  plotting  routine 
for  the  Cyber  750  machine.  A  mean  displacement  curve  was  then  obtained 
for  each  set  of  creep  tests  by  averaging  the  displacement  vs.  time  curve 
of  each  test  phase.  Next,  all  specimens  were  dissected  and  measured  to 
obtain  data  needed  in  formulating  a  finite  element  model.  Iri  the  final 
phase  of  the  analysis,  a  viscoelastic  ax i symmetric  fwiite  element,  model 
was  used  to  quantify  the  exper  imen  ta  1  data.  A  three  parameter  Kelvin 
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solid  was  employed  in  the  finite  element  solution  method./  The  result  of 
this  paper  will  be  used  to  construct  a  dynamic  model  of  the  vertebral 
column . 


CHAPTER  1 

INTRODUCTION 

1.1.  Background 

Numerous  studies  of  the  intervertebral  disc  of  both  an  experimental 
and  an  analytical  nature  have  been  done.  Although  the  work  presented 
here  was  done  on  the  rhesus  monkey  (Macaca  mulatta),  numerous  comparisons 
and  correlations  can  be  made  between  it  and  those  made  on  human 
specimens.  Work  done  by  Kazarian  (8)  addresses  the  specific  subject  of 
comparisons  between  different  species. 

Several  experimental  papers  have  been  concerned  with  the  fact  that 
under  some  loading  conditions  the  intervertebral  disc  behaves  in  a  visco¬ 
elastic  manner.  Several  of  the  early  papers  studying  the  time  dependent 
observation  can  be  attributed  to  Virgin  (17)  3nd  Hirsch  (7).  Markolf 
(15)  considered  several  cyclic  loads  on  discs  in  varying  states  of 
surgical  herniation  and  always  found  the  disc  to  behave  viscoel astically. 
More  recently,  people  have  been  concerned  not  only  with  experimentation 
but  also  with  analysis  of  the  data  either  by  modeling  of  the  results  or 
by  altering  the  specimen  and  roperforming  the  experiment.  Lin  and  Liu 


(14)  considered  several  types  of  loading  patterns  and  concluded  that  during 
pure  axial  loading  the  posterior  elements  play  a  minor  role  in  load  bearing. 
Nachemson  (16)  felt  that  a  gross'y  herniated  disc  did  not  exhibit  a 
significantly  different  response  from  the  healthy  specimens.  Kazarian  (11), 
who  did  numerous  compressive  tests,  concluded  that  under  certain  conditions 
the  disc  would  not  behave  in  a  viscoelastic  manner. 

Although  the  experimental  studies  showed  viscoelastic  phenomena,  the 
analytical  methods  approach  modeling  from  a  variety  of  viewpoints. 
Belytschko  did  an  axisymmetric  finite  element  analysis  using  an  ortho¬ 
tropic  annulus  fibrosus  with  an  incompressible,  inviscid  nucleus  pulposus. 
Kulak  (12)  has  also  done  nonlinear  modeling  using  finite  elements.  He  did 
his  work  based  on  analysis  by  Fung  (5)  which  stated  that  stress  in  some 
biological  materials  can  be  broken  down  into  an  elastic  and  history 
dependent  part.  Lin  (13)  did  a  three  dimensional  finite  element  model 
assuming  orthotropic  material  properties  for  the  disc  and  vertebrae.  Burns 
(3)  discusses  the  idea  of  viscoelastic  modeling  using  Kelvin  solid  units  and 
came  to  the  conclusion  that,  the  three-unit  solid  results  in  the  best 
quantitative  results.  Hinrichsen  (6)  employed  an  axisymmetric  viscoelastic 
finite  elemient  model  using  a  three  parameter  Kelvin  solid. 

1.2  Purpose 

In  aircraft  flight,  the  body  can  be  subjected  to  severe  physical 
loads;  emergency  egress,  high  G  turns,  and  random  vibrations.  The 
vertebral  column  must  absorb  much  of  this  loading.  As  a  result, 
intervertebral  disc  herniation  sometimes  occurs.  In  this  light,  this 
thesis  will  investigate  the  influence  of  surgical  herniation  on  the 
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viscoelastic  properties  of  the  intervertebral  disc.  Both  experimental  and 
analytical  phases  will  be  carried  out  in  this  study. 

1.3  Anatomy 

A  brief  discussion  of  rhesus  monkey  spinal  anatomy  will  be  helpful  in 
understanding  the  remainder  of  this  thesis.  First,  the  gross  anatomy  of 
the  spine  will  be  examined  followed  by  a  more  specific  discussion  of  the 
vertebral  segment. 

The  vertebral  column  is  divided  into  thirty-three  vertebrae:  7 
cervical,  12  thoracic,  7  lumbar,  6  sacrum,  and  4  coccyx  (see  fig.  (1.1)). 
In  addition,  most  vertebrae  have  a  disc  between  two  adjacent  vertebrae 
called  the  i ntervertebral  disc.  Several  important  aspects  of  the 
vertebral  segment  are:  posterior  elements,  longitudinal  ligaments, 
vertebral  centrum  and  disc  (see  fig.  (1.2)). 

The  concern  of  this  work,  however,  is  mainly  in  characteristics  of 
the  vertebral  segment  which  is  made  ip  of  two  parts:  vertebral  centrum 
and  intervertebral  joint. 

The  vertebral  centrum  is  composed  of  three  types  of  hard  tissue: 
cortical  bone,  trabecular  bone,  and  bony  endplate  (see  fig.  (1.4)).  The 
cortical  bone  is  a  very  thin  and  strong  material  which  surrounds  the  outer 
surface  of  the  centrum.  The  trabecular  bone  is  a  soft  material  with  a 
matrix  very  similar  to  that  of  sponge.  The  bony  endplate  is  a  hyaline 
cartilage  structure  that  separates  the  disc  from  the  centrum. 

The  disc  also  has  several  components:  nucleus  pulposus  and  annulus 
fibrosis  (sea  fig.  (1.3)).  The  nucleus  makes  up  approximately  half  of  the 
planar  cross-section  of  the  disc  an-1  is  composed  of  fibrous  strands  in 
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fig.  1.2  (from  13)  Vie*  of  vertebral  unit  with  posterior  elements 
attached. 


1.3  Agood  illustration 
the  disc 


of  the  complexity  of 


region 


a  muco -protein  gel  with  a  water  content  of  70-90%.  The  annulus  is  a 
relatively  complicated  structure.  It  consists  of  fibrous  tissue  in 
concentric  laminate  bands.  All  fibers  in  any  one  band  have  the  same 
orientation  but  opposite  to  that  of  its  two  adjacent  bands.  An  analogous 
structure  to  the  disc  is  the  jelly  donut  where  the  nucleus  is  the  jelly 
and  the  annulus  serves  as  the  dough. 

1.4  Assumptions  and  General  Approach 

The  general  approach  of  this  thesis  can  be  divided  into  two  sections: 
experimentation  and  analysis. 

The  experimentation  involved  tests  of  lumbar  intervertebral  disc  at 
the  Lj-U  level  in  varying  states  of  surgical  herniation.  At  all 
stages  of  herniation,  a  creep  test  was  done  with  an  eight-hour  compressive 
load  followed  by  a  sixteen -hour  relaxation  time.  The  experimental  results 
indicated  that  the  creep  behavior  of  the  specimen  is  not  altered  unless 
the  disc  is  grossly  herniated. 

The  analysis  of  the  experimental  data  was  done  in  the  following  manner. 
The  specimen  was  modeled  assuming  rotatioria1  symmetry  about  the  center  of 
the  disc.  A  finite  element  analysis  was  done  to  obtain  the  material 
properties  of  the  disc  when  the  disc  is  modeled  as  a  spring  in  series  with 
a  Kelvin  solid.  In  addition,  stress  profiles  ana  deformed  plots  were 

obtained  of  the  specimen  at  various  times  and  stages  of  herniation. 

The  final  conclusion  indicated  that,  in  all  stages  of  surgical  hernia¬ 
tion,  the  disc  behaves  in  a  viscoelastic  manner.  Moreover,  the  larger 
degree  of  herniation  cor  responded  wel 1  with  a  larger  initial  displacement  and 
a  faster  creep  rate  over  lime. 
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CHAPTER  2 

EXPERIMENTATION 


2.1  Introduction 

The  primary  aspects  of  the  experimentation  were:  considerations, 
protocol  and  results.  In  the  considerations,  the  experimental  data  must, 
conform  to  an  established  method  of  analysis.  In  addition,  the  experi¬ 
mental  protocol  must  be  structured  so  that  all  test  specimens  are  tested 
as  equally  as  possible.  A  thorough  experimental  procedure  is  given  in 
Appendix  A. 

2.2  Considerations 

In  determining  what  types  of  information  were  desired  from  experi¬ 
mentation,  several  factors  had  to  be  considered.  These  considerations 
included  shape,  movement  during  testing,  age  and  condition,  drying  out, 
and  method  of  herniation  of  the  specimen. 

Several  factors  were  important  in  determining  shape  of  the  specimen. 
It  should  have  symmetry  with  respect  to  its  top  and  bottom  half.  In 
addition,  the  shape  of  the  specimen  must  lend  itself  to  an  axi symnetric 
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analysis.  These  two  considerations  were  sat isified  by  using  the  disc  at 


the  L-j -1~2  level  of  the  rhesus  monkey  spine. 

Another  cons iderat ion  which  the  experimentation  must  meet  is  the 
movement  of  the  specimen  must  be  in  a  purely  compressive  fusion  to  insure 
that  an  axisymmetric  analysis  is  justified.  Two  pretests  were  run  on 
sample  specimens  to  determine  if  pitch  or  torsion  of  the  top  centrum 
relative  to  the  bottom  centrum  occurred  (see  fig.  (2.1)  for  examples  of 
pitch  and  torsion).  Both  pretest  specimens  displayed  pitch  and  torsion. 
To  insure  against  the  non-vertical  movement  during  the  experimentation, 
each  specimen  was  cupped  with  a  dental  acrylic.  The  cupping  also  allowed 
for  a  more  evenly  distributed  loaa  across  the  top  edge  of  the  test 
sample.  Finally,  to  measure  any  pitch  or  torsion  which  might  have 
occurred,  each  specimen  had  two  pins  placed  in  each  centrum  and  time 
delay  photographs  were  taken  to  measure  the  pitch  and  torsion  (see  fig. 
(2.2  and  2.3)). 

Age  and  condition  of  the  specimen  were  also  an  important  factor 
before  the  experimentation.  If  either  of  these  two  varied  over  a  wide 
range,  difficulty  could  arise  in  the  analysis  because  both  age  and 
condition  affect  the  strength  characteristics  of  the  vertebral  column. 
In  addition,  one  could  have  difficulty  in  obtaining  a  mean  when  these  two 
factors  run  the  gamut.  To  determine  age  and  condition  of  each  specimen, 
each  vertebral  column  was  X-rayed  in  the  anterior-posterior  and  the 
sagittal  planes.  The  age  of  ttie  specimen  was  then  determined  by  a 
veterinarian  by  measuring  the  thickness  of  the  epiphyseal  plate.  All 
spines  were  determined  to  have  come  from  young  adult  rhesus  monkeys.  One 
specimen  was  found  to  be  diseased  and  as  a  result  the  data  gathered  from 
it  was  not  used  in  obtaining  mean  displacement  curves. 


fig  2.2  "A"  is  a  picture  of  the  uncupped  prepared  specimen.  "B 
is  a  view  of  the  cupped  specimen  ready  for  testing.  Note  that 
both  specimens  have  two  pins  in  each  vertebral  centrum.  Ihese 
pins  are  used  to  measure  pitch  and  torsion  from  the  time  delay 
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Since  the  testing  was  carried  out  over  a  period  of  twenty-four 
hours,  drying  out  of  the  specimen  was  considered.  In  addition,  one  is 
concerned  with  the  specimen  absorbing  too  much  fluid  because  altering 
water  content  of  the  disc  has  been  shown  by  Anderson  (1)  to  change 
material  properties.  To  insure  against  either  of  these,  several  steps 
were  taken.  First,  during  the  experimentation,  the  specimen  was  placed 
in  an  environment  of  approximately  100%  relative  humidity.  And  during 
the  storage,  the  specimens  were  placed  in  a  saline  rich  environment  and 
kept  refrigerated. 

The  final  concern  was  determining  the  method  of  herniation.  Ideally 
one  would  like  to  remove  the  nucleus  pulposus  from  the  disc  without 
altering  the  annulus  fibrosis.  In  the  first  attempt  at  herniation,  a 
hypodermic  needle  is  used  to  remove  the  nucleus  pulposus.  The  experi¬ 
mental  results,  however,  indicated  difference  from  the  healthy 
specimens.  In  a  second  c.  1 1  e  t  p :  te  herniate  a  t.reaph  ine,  a  surgical 
instrument  similar  in  stnict.u--  t<  a  t >  ■ 1 .  s  *,  was  used  to  remove  a  small 
piece  of  the  annulus  and  thus  neiia,  .  leus  as  a  structural 

component  of  the  disc.  At  thit  p  < .  i  ■,  t  .  it  was  necessary  to  make  the 
assumption  that  as  the  nucleus  passed  through  the  small  hole  in  the 
annulus,  made  by  the  trephine,  it  contributed  to  the  load  resistance  in 
the  annulus  and  thus  created  no  discontinuity. 

2.3  Protocol 

The  exper imental  protocol  was  divided  into  three  parts: 

Preparation,  carrying  out  of  load  cycles  in  each  stage  of  herniation,  and 
dissection  of  the  specimen. 
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Preparation  of  the  specimen  involved  several  steps.  The  specimen 


was  first  removed  from  the  L^-L^  level  of  the  vertebral  column  of 
young  adult  rhesus  monkeys.  Next,  all  soft  tissue  was  removed  from  both 
the  centrums  and  the  posterior  elements.  Care  was  taken  not  to  damage 
the  disc  and  the  anterior  and  the  posterior  longitudinal  ligaments. 
The  ligaments  were  kept  intact  because  to  remove  them  would  have  resulted 
in  damage  to  the  disc.  The  final  step  before  cupping  of  the  specimen  was 
to  remove  the  posterior  elements.  The  cupping  of  each  specimen  was  done 
with  a  minimal  amount  of  dental  acrylic  to  insure  the  cup  was  not  clamped 
to  the  specimen  and  thus  possibly  create  a  moment  during  the  loadings. 
In  addition,  a  special  jig  was  used  to  insure  both  cups  were  set  parallel 
to  one  another. 

The  first  load  cycle  involved  testing  of  the  healthy  specimens.  The 
tests  were  carried  out  in  specially  built  humidity  chambers  to  insure  the 
water  content  of  the  specimen  did  not  change  (see  fig.  2.3).  The  dis- 
placment  of  the  top  edge  of  the  specimen  was  recorded  through  all  load- 
unload  cycles  with  an  LVOT  (Linear  Variable  Differential  Transducer) 
whose  accuracy  was  0.0002  centimeters.  In  addition,  the  LVDT  was 
calibrated  and  set  to  zero  before  each  loading.  Throughout  the  test, 
time  delay  photographs  of  the  specimens  were  taken  every  twenty-four 
seconds  to  obtain  data  on  pitch,  torsion,  and  bulge  data  of  the  disc. 
After  this  preparation,  the  specimens  were  then  subjected  to  an 
eight-hour  load  phase  and  sixteen-hour  unload  phase  in  the  test  cycle. 
During  the  load  phase,  a  total  dead  weight  of  6.7  kiloponds  was  applied. 
During  the  unload,  the  weight  was  removed  to  obtain  the  relaxation  curve. 

In  the  second  load  cycle,  the  general  procedure  was  the  same  as  the 
first  except  an  attempt  was  made  to  herniate  the  disc  with  fourteen- gauge 


hypodermic  needle.  The  nucleus  was  then  drawn  out  using  suction. 

In  the  final  load  cycle,  a  trephine  was  used  to  herniate  the  disc. 
In  this  method,  the  trephine  made  a  5mm  hole  in  the  annulus  allowing  the 
nucleus  to  flow  out  through  the  annulus  during  the  loading. 

After  the  fina1  loading  test,  all  specimens  were  dissected  to  obtain 
measurements  needed  in  formulating  the  finite  element  model. 

2.4  Results 

Not  only  did  the  experimental  results  conform  with  expectations,  but 
the  graphed  data  gave  a  clear  indication  of  the  method  of  analysis.  In 
examination  of  the  time  delay  photographs,  all  specimens  showed  a  purely 
compressive  displacement  with  no  measurable  pitch  or  torsion.  In 
addition,  the  disc  was  shown  to  bulge  two  to  four  percent  during  the  load 
cycles.  The  creep  curves  of  displacement  versus  time  in  the  following 
pages  were  established  by  obtaining  the  mean  of  all  five  specimens  over 
time.  The  last  curve,  which  shall  be  referred  to  as  osteo,  is  that  of  a 
diseased  specimen  and  will  be  examined  separately  from  the  other  data. 

The  osteo  specimen  was  diagnosed  by  a  veterinarian  to  have 
osteochondritis.  From  a  physical  standpoint,  this  disease  may  be 
explained  by  the  fact  that  the  bony  endplate  extended  into  the  nucleus 
pulposus  to  the  point  where  the  disc  had  almost  no  nucleus  pulposus.  The 
value  of  this  specimen  lies  in  comparing  its  material  properties  to  that 
of  the  healthy. 

The  mean  displacement  versus  time  curve  for  the  first  two  test 
cycles  turned  out  to  be  identical  (see  fig.  (2.4)).  Several  factors  con- 


tribute  to  this  similarity.  When  the  disc  was  herniated  with  the  needle, 
almost  no  fluid  was  extracted  from  nucleus.  As  a  result,  one  would 
assume  very  little  or  no  change  in  the  creep  curve.  A  second  reason  for 
the  lack  of  change  between  the  two  curves  was  a  self -sealing  mechanism 
proposed  by  White  (18)  that  states  the  annulus  will  reseal  itself  to 
prevent  escape  of  the  nucleus  when  the  disc  is  damaged.  Because  the 
curve  for  the  first  two  load  cycles  was  identical,  the  results  from  the 
first  herniation  shall  not  be  further  analyzed.  Finally,  the  mean 
displacement  versus  time  curve  shows  the  specimens  behaved  in  a  visco¬ 
elastic  manner. 

The  second  experimental  curve  (see  fig.  (2.5))  which  is  the  mean 
curve  for  those  specimens  herniated  with  a  trephine  shows  marked 
differences  from  the  previous  curve.  Not  only  is  the  initial  displace¬ 
ment  larger  than  the  healthy  specimen,  but  the  grossly  herniated  one- 
shows  a  faster  creep  rate  over  time.  The  grossly  herniated  specimen  also 
displays  viscoelastic  characteristics  over  time. 

The  final  curve  of  the  osteo  specimen  (see  fig.  (2.6))  is  perhaps 
the  most  confusing.  The  creep  curve  shov/s  a  smaller  initial  displacement 
than  the  healthy  specimen.  In  addition,  the  creep  rate  for  this  specimen 
was  faster  than  that  of  the  healthy.  Once  again,  the  displacement  versus 
time  for  this  curve  exhibits  creep. 

The  preliminary  results  to  this  point  are  excellent.  Even  though 
the  first  attempt  at  herniation  was  not  successful,  the  results  were 
easily  explained  and  understood.  Secondly,  the  curves  obtained  from  this 
experiment  are  very  similar  to  those  obtained  by  Virgin  (17),  Markolf 
(19),  and  Nachemscn  (16).  In  the  next  chapter,  a  viscoelastic  theory 
will  be  developed  to  assist  in  modeling  of  the  data. 


16 


.GO  35. 


Experimental  results  of  fully  herniated. 
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The  two  units  may  be  combined  either  in  series  as  a  Maxwell  model 


or  in  parallel  as  a  Kelvin  solid. 
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The  Maxwell  model  forms  on  the  addition  of  two  strains. 
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The  constitutive  equation  for  the  Kelvin  solid  is  developed  in  an 
analogous  manner  based  on  the  sum  of  two  stresses 
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which  is  now  an  expression  for  the  creep  strain 
implemented  into  the  following  form 

let  n-=%,  E-  6-/ 


r 


d  ^  ~  dx 

%  i  L 


3.3  Three-Dimensional  Viscoelastic  Theory 
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Three-dimensional  viscoelastic  theory  is  based  on  two  ideas:  the 
constitutive  equation  of  one  dimensional  N  Kelvin  units  model  and  the 
viscoelastic  law.  In  a  stress  analysis  of  a  three-dimensional  element, 
the  stress  tensor 
^~xy  f'yz 
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can  be  broken  down  into  two  components:  A  trace  matrix  and  a  deviatorvic 
matrix.  In  a  similar  manner,  the  strain  tensor  can  be  broken  into  two 
matrices:  dilatation  matrix  which  is  volumetric  change  without  change  of 
shape,  and  a  distortion  matrix  which  is  change  of  shape  at  constant 
volume. 
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The  first  matrix  on  the  right  hand  side  is  the  dilatation  and  the  second 
is  the  distortion  matrix. 

The  viscoelastic  law  uses  the  constitutive  equation  for  the  Kelvin 
solid  model  and  the  idea  of  separating  the  stress  and  the  strain  matrices 
into  two  components  to  form  a  three  dimensional  viscoelastic  equation. 
If  one  has  an  isotropic  viscoelastic  material,  the  hydrostatic  stress,  s, 
must  correspond  to  a  dilatation,  e,  with  no  distortion.  Using  this  idea. 


the  following  equation  is  obtained 
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The  second  equation  from  the  viscoelastic  law  yields  a  relation  between 
the  deviatorvic  stress  tensor  and  the  distortion  matrix 
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where  S  &  E  correspond  to  the  components  of  the  stress  and  strain  devi¬ 
ators.  The  four  operators  P",  Q",  P'  arid  Q'  are  entirely  independent  of 
one  another.  In  order  to  find  solution  to  the  above  equations  given  a 
body  similar  to  the  intervertebral  disc  would  be  extremely  difficult  if 
at  all  possible.  For  this  reason,  one  is  driven  to  the  finite  element 
method. 
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3.4  Elastic  Analysis 


The  elastic  analysis  of  an  axisymmetric  problem  is  divided  into 
several  parts.  The  strain  is  expressed  by  (see  fig.  3.1  for  coordinates) 
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or  in  term  of  the  nodal  displacment 
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and  A  is  the  area  of  the  element, 
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Similarly  the  stress  is  expressed  by 
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where  [D]  is  the  elastic  material  property  matrix 
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The  overall  solution  equation  to  this  problem  is 

[kJH  =  [p]  r  fei  Ci  -zs) 

where  [K]  is  the  elastic  stiffness  matrix,  {p}  is  the  applied  load 
vector,  and  ^q]  is  the  applied  load  due  to  initial  strains.  In  integral 
form,  [K]  and  may  be  represented  as 
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The  solution  method  for  this  problem  was  formulated  assuming  a  centroidal 
radius,  r,  and  centroidal  height,  z,  where 
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The  result  is  a  new  [R]  matrix  which  shall  be  referred  to  as  TB]. 
At  this  point  one  can  numerically  obtain  the  stiffness  matrix  [K]  and  the 
residual  strain  matrix  Q  without  integration.  The  new  expressions  for 
[Kl  and  Q  are 
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At  this  point,  one  can  obtain  the  nodal  displacements  through  the 
solution  of  equation  (3-7). 

3.5  Viscoelastic  Analysis 


The  viscoelastic  analysis  involves  the  combined  concept  of  elastic 
strain  in  addition  to  the  incompressibility  of  the  time  dependent 
char acter i st i cs  .  Thus,  the  viscoelastic  strain  can  be  determined  by 
recognizing  that  the  incorporation  of  the  distortion  strain  tensor 
assumes  zero  volummetric  change.  The  change  in  volume  is  contained  in 
the  elastic  relationships.  In  order  to  evaluate  the  two-dimensional 
viscoelastic  stress,  the  strain  is  incrementally  determined  in  the  same 
manner  shown  in  equation  (3-13).  for  the  finite  element  formulation,  the 
strain  increment  is  expressed  by 
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In  order  to  use  each  strain  increment the  residual  strain  matrix  Q 
is  altered  for  each  time  stepYi'--.  As  a  result,  one  can  formulate  a 
viscoelastic  problem  by  using  incremental  time  steps  and  strain  rates. 
The  flow  chart  for  the  program  used  in  the  analysis  is  given  in  fig.  3.2. 
Thus  it  is  possible  with  the  above  equation  to  incorporate  the  three 


parameter  model  into  the  analysis.  Therefore,  the  strain  tensor  can  be 


written  as 
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where  C  is  the  elastic  strain  tensor  which  includes  volummetric  and 

ve 

distortion  deformation  and  e.  is  the  viscoelastic  component  for  the 
incompressible  strain  present  in  equations  (3-35)  and  (3-36)  where 


CHAPTER  4 


FINITE  ELEMENT  ANALYSIS  OF  EXPERIMENTAL  DATA 

4.1  Introduction 

The  analysis  of  the  experimental  data  is  divided  into  three  parts: 
formulation  of  a  finite  element  model,  coarse  mesh  and  fine  mesh 
analysis.  The  three  parameter  Kelvin  solid  will  be  used  to  model  the 
experimental  data  (see  fig  4.1).  The  three  material  constants  Eq, 
qQ,  and  q^  may  be  interpreted  in  the  following  manner.  The  elastic 
spring  constant,  Eq  is  used  to  match  the  initial  displacement  of  the 
specimen.  Whereas,  the  second  spring  constant,  qQ,  determines  the 
displacement  over  time.  The  final  constant,  q^,  influences  the  slope 
of  the  creep  curve. 

4.2  Formulation 

The  formulation  of  model  has  three  key  sections:  reasons  for  an 
axisymmetric  model,  modeling  of  the  vertebral  centrum  and  finally 
modeling  of  the  intervertebral  disc. 

The  axisymmetric  analysis  is  an  excellent  form  of  analysis,  when 
studying  the  intervertebral  disc,  for  several  reasons.  In  a  top  view  of 
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fig.  4.1  A  three  parameter  kelvin  so}  id. 


fig.  4.2  The  dark  line  illustrates  the  shape  of  the  disc  from  the 
top.  The  dotted  line  is  the  equivalent  circular  representation. 
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the  specimen  (see  fig.  (4. "I)),  the  specimen  itself  is  very  close  to  a 
circle  in  shape.  In  addition,  the  axisymmetric  model  has  been  used  in 
the  past  by  several  different  investigators;  Belytschko  (2),  Kulak  (12), 
and  Hinrichsen  (6).  Finally,  this  type  of  analysis  enables  one  to  obtain 
an  understanding  of  the  stresses  in  the  disc  without  the  problems  of  a 
three-dimensional  analysis.  For  the  same  level  of  depth  of  analysis  in 
terms  of  computer  time  and  core  memory,  one  can  model  the  interface 
between  the  disc  and  the  centrum  much  more  accurately  in  an  axisymmetric 
than  in  a  three-dimensional  analysis.  Now  that  the  form  of  analysis  has 
been  obtained,  the  next  step  is  to  examine  the  centrum  and  the  disc 
separately. 

Two  problems  existed  in  modeling  of  the  vertebral  centrum:  Tirst, 
determining  material  properties  for  trabecular  bone,  cortical  bone,  and 
bony  endplate,  and  secondly,  accurately  modeling  the  geometry  of  the 
bone.  The  modulus  of  elasticity  and  Poisson's  ratio  for  human  cortical 
and  trabecular  bone  can  be  obtained  from  Belytschko  (2).  In  addition, 

the  modulus  of  elasticity  of  the  vertebral  centrum  for  both  the  human  and 

rhesus  monkey  are  known  from  Kazarian  (8,9).  However,  the  material 
constants  for  cortical  or  trabecular  bene  of  the  rhesus  are  not  known. 
Kazarian  (8)  discusses  the  comparing  of  material  characteristics  between 
species.  To  obtain  the  material  constants  needed  for  the  finite  element 
model,  several  assumptions  were  made.  First,  it  was  assumed  Poisson's 
ratio  for  human  and  rhesus  bones  were  the  same.  Secondly,  since  centrum 
properties  were  available  from  both  species,  the  ratio  of  the  moduli  of 

elasticity  was  used  to  estimate  a  modulus  for  the  cortical  and  a  modulus 

for  the  trabecular  bone  of  the  rhesus  monkey.  The  ratio  used  to 
determine  these  constants  was 
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The  final  material  properties  were:  cortical  =  6.207x10  Kp/cm  , 

F  T  2 

trabecular  =  2.80x10  Kp/cm  and  both  used  \j  =  0.25.  The  bony 
endplate  material  properties  are  not  clearly  known.  However,  it  is  known 
to  be  harder  than  trabecular  bone.  As  a  result,  the  bony  endplate  was 
modeled  with  the  same  material  properties  as  the  cortical  bone.  The 
shape  of  the  vertebral  centrum  also  presented  some  problems.  To  solve 
th^s  difficulty  when  the  specimens  were  dissected,  all  measurements 
necessary  to  formulate  the  finite  element  model  were  taken.  It  was  found 
that  both  the  bony  endplate  and  the  cortical  bone  were  extremely  thin. 
As  a  result,  a  very  fine  mesh  is  needed  to  model  both  sections 
adequately.  In  addition,  the  fine  mesh  is  needed  to  provide  a  good 
transition  of  the  load  from  the  trabecular  bone  to  the  disc.  The  disc 
will  be  modeled  as  a  homogeneous  body  for  several  reasons.  First,  since 
the  annulus  is  an  orthotropic  material  and  since  it  clearly  displayed 
viscoelastic  properties  in  the  final  creep  test,  to  model  it  thoroughly 
would  involve  an  orthotropic  viscoelastic  analysis  with  nine  material 
constants.  Obviously  determining  nine  independent  material  constants 
is  an  awesome  task.  A  second  problem,  in  dividing  up  the  disc,  is  how 
does  one  adequately  model  the  fluid  movement  in  the  nucleus,  annulus,  and 
between  the  two  materials?  The  problem  here  again  is  very  difficult; 
fluid  movement  in  a  deforming  body.  In  modeling  the  disk  as  a  unit,  one 
avoids  many  of  the  problems  mentioned  above.  However,  one  does  obtain 
information  of  the  stress  distributions  and  displacement  profiles. 
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4.3  Coarse  Mesh  Analysis 


The  coarse  mesh  analysis  served  to  provide  excellent  initial  guesses 
to  the  fine  mesh.  However,  from  a  finite  element  standpoint,  the  coarse 
mesh  had  several  short  comings. 

The  primary  use  of  the  coarse  mesh  was  to  provide  information  needed 
for  the  fine  mesh.  Since  the  coarse  mesh  had  approximately  one-fifth  the 
number  of  variables  of  the  fine  mesh,  one  could  run  many  more  coarse  mesh 
jobs  in  a  day  than  fine  mesh  jobs.  In  addition,  a  clear  relation  existed 
between  the  spring  constants  (EQ,qo)  for  the  coarse  model  and  the 
constants  for  the  fine  mesh  (see  fig.  (4.3)).  As  a  result,  one  could 
narrow  in  on  material  constant  for  the  fine  mesh  quickly.  A  second 
advantage  to  the  coarse  mesh  is  that  it  provided  an  almost  exact  value 
for  the  dashpot  constant,  for  the  fine  mesh. 

However,  the  coarse  mesh  did  have  several  drawbacks.  The  most 
evident  one  is  the  high  aspect  ratios  in  both  the  cortical  bone  and  bony 
endplate  region  (see  fig.  (4.3)).  The  second  major  problem  to  this  mesh 
was  the  lack  of  transition  from  the  trabecular  to  the  disc  region.  The 
utility  of  the  coarse  mesh  was  in  its  ability  to  provide  a  good  initial 
guess  for  the  fine  rresh. 

4.4  Fine  Mesh  Analysis 

The  fine  mesh  analysis  will  be  broken  down  into  several  sections. 
First  a  discussion  of  the  formulation  of  the  fine  mesh,  next  each  set  of 
experimental  data  will  be  analyzed  individually.  First  the  healthy,  then 
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fig.  4.4  The  above  graph  illustrates  the  relation  between 
the  spring  constants,  and  q^,  for  the  coarse  to  the 
fine  mesh. 
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the  fully  herniated,  and  finally  the  osteo  specimen  will 
In  the  last  section  of  this  chapter,  a  comparison  will  be  made  between 
the  three  analyses. 

The  fine  mesh  model  was  obtained  from  two  parameters:  the  need  to 
show  transition  from  one  material  to  another  and  secondly  using  a 
progressively  finer  mesh  to  obtain  a  converged  solution.  The  transition 
from  one  material  to  another  was  the  primary  reason  for  the  fine  mesh. 
The  aspect  ratio  for  all  of  the  elements  was  reduced  until  the  maximum 
was  3  to  1.  At  this  point,  the  mesh  was  fine  enou^i  to  obtain  a 
converged  solution  also. 


i 
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4.4.1  Analysis  of  the  Healthy  Specimen 

The  analysis  of  the  healthy  data  can  be  examined  at  two 
levels:  the  ability  of  the  model  to  match  the  displacment  of  the  specimen 
and  secondly,  a  disucssion  of  stress  flow  in  the  specimen  during  the 
creep  test. 

The  finite  element  model  gave  an  excellent  fit  to  the 
experimental  results.  For  the  fine  mesh,  the  material  constants  which 
gave  the  best  fit  were  Eo=  36.50  *(D/cm^,  qQ  =16.21  ^/cm^, 

=44060  kP/cm  sec,  and  =0.48.  The  actual  curve  fit  is  shown  in 
fig.  4.6.  In  addition,  the  radial  displacment  of  the  disc  was  very  close 
to  the  experimental  as  can  be  seen  in  fig.  4.7.  One  can  also  see  how  the 
model  deformed  over  time  by  examining  figs.  4.8-4.10.  Please  note  the 
displacements  in  the  figures  are  multiplied  by  a  factor  of  two  to  more 
clearly  illustrate  the  deformations. 

The  finite  element  model  also  gave  a  good  indication  of  the 
stress  within  the  specimen.  The  graph  of  the  0,  hoop  stress  versus 
radius  in  the  centrum  (fig.  4.11)  shows  that  the  cortical  bone  is  acting 
like  shell  in  compression  restraining  the  trabecular  bone.  This  becomes 
obvious  if  one  looks  at  the  geometry  of  the  cortical  bone.  It  is  in  the 
shape  of  an  inwardly  curved  surface  referred  to  as  an  anticlastic 
surface.  Thus  the  radial  tension  and  axial  compression  produce 
compression  in  the  hoop  direction.  One  can  also  note  that  the  stress 
increases  with  time.  Since  the  centrum  has  the  same  material  properties 
in  all  three,  analyses,  it  will  not  be  covered  in  the  discussion  of  the 
fully  herniated  or  the  osteo  analysis.  The  stress  data  also  gave 
numerous  indications  of  the  stress  flow  in  the  disc.  The  graph  of  the 


radial  stress  versus  radius  (Fig.  4.12)  clearly  shows  the  hony  endplate 
is  compressing  the  disc  except  in  the  outer  area  of  disc  where  the  bony 
endplate  does  not  cover  the  disc.  In  addition  one  can  see  that  the 
stress  level  increases  with  time.  The  second  graph  in  the  disc  region  of 
the  vertical  stress  versus  radius  (fig.  4.13)  again  shows  the  bony 
endplate  is  placing  the  disc  in  compression.  One  can  again  see  the 
stress  increases  with  time.  The  final  pertains  to  the  hoop  stress  in  the 
disc  (fig.  4.14)  and  shows  an  increase  in  the  stress  level  with  time. 
The  graph  also  illustrates  that  the  outer  disc,  which  is  not  covered  by 
the  bony  endplate,  is  in  tension.  This  indicates  the  disc  annulus  region 
is  resisting  load  as  a  shell  structure.  The  outward  radial  curvature  is 
indicative  of  a  tension  stress  in  the  hoop  direction  arising  from  the 
axial  compressive  stress. 

4.4.2  Analysis  of  the  Fully  Herniated  Specimen 

The  analysis  of  the  fully  herniated  test  results  will  be 
covered  in  the  same  manner  as  the  healthy.  The  displacements  will  be 
discussed  first  followed  by  the  stress  analysis. 

The  finite  element  model  once  again  gave  an  excellent  match 
to  the  experimental  curves.  The  material  constants  which  gave  the  best 
fit  were  Eq-  20.61  Kp/cm^,  qQ^  5.97  Kp/crr^,  q^=  18500  ^P/crn  sec, 
and  =  .480.  A  graphical  illustration  of  the  fit  is  given  in  fig.  4.15. 
One  can  also  see  deformation  of  the  model  in  figs.  4.16  -  4.17. 

The  stress  profiles  also  give  a  good  indication  of  the  stress 
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fig.  4.C  Finite  element,  solution  for  vertical  displacement 
of  tile  healthy  specimen 


.29  14Z.B6 


fig.  A . 8  fine  finite  model  of  healthy 
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fig.  4.12  Radial  sire sr.  versus  radius  in  the  disc  of  the 
healthy  specimen. 
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redistribution  with  time.  The  first  graph  in  fig.  4.18  of  the  radius 
stress  in  the  disc  shows  the  stress  decreasing  as  the  radius  increases. 
In  addition,  one  can  see  the  total  radial  stress  increases  with  time. 
The  second  graph  in  the  disk  region  of  the  vertical  stress  versus  radius 
(fig.  4.19)  again  shows  the  bony  endplate  is  placing  the  disc  in 
compression.  However,  in  this  curve  the  entire  disc  is  in  compression. 
One  can  again  see  the  stress  increases  with  time.  The  final  curve 
pertains  to  the  hoop  stress  in  the  disc  (fig.  4.17)  and  shows  an  increase 
in  stress  level  with  time.  This  hoop  stress  profile  is  again  showing 
that  the  nature  of  restraint  by  the  disc  as  a  shell  structural  element. 

4.4.3  Analysis  of  the  Osteo  Specimen 

Before  discussing  the  analytical  results,  a  brief  description 
of  the  diseased  state  of  the  osteo  specimen  will  be  helpful.  The  osteo 
specimen  was  diagnosed  by  a  veterinarian  to  have  osteochondritis.  From 
the  X-ray,  one  could  see  the  bony  endplate  of  the  upper  centrum  extended 
well  into  the  nucleus  to  point  at  which  almost  no  nucleus  was  present  in 
the  disc.  The  same  analytical  model  was  used  to  study  this  specimen  as 
the  previous  two  data  sets.  In  the  analysis  of  the  osteo  experimental 
data,  the  same  type  of  information  was  obtained  as  in  the  two  previous 
models.  The  displacement  curves  will  be  studied  first,  followed  by  the 
stress  analysis. 

The  finite  element  model  once  again  gave  a  very  close 
correlation  to  the  experimental  data.  The  best  curve  fit  were  with 
material  constants  of  Eq=  56. 25  Kp/cm2  qQ=  9.218  Kp/cm?  q^- 
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fig.  4.16  Displacement  of  fully  herniated  specimen  at  t-=0 


RflOi 


fig.  4.19  Vertical  stress  versu> radius  jn  the  disc  of 
the  fully  herniated  specimen. 
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f'iy*  4.20  Hoop  stress  versus  radius  in  trie  disc  of  the 
fully  herniated  specimen. 


*  u  »iwi  *r 


0G  a. 1  1  0.23  0.34  0.46  0.57  0.69  0.80 

RADIUS  (CM) 


27030  ^/cm  sec  and  =  .480.  The  finite  element  solution  is  shown  in 
fig.  4.21.  Once  again  deformed  plots  were  made  of  the  specimen  and  are 
shown  in  fig.  4.22-4.23. 

The  stress  profiles  were  made  to  show  the  stress  redistri¬ 
bution.  The  radial  stress  (fig.  4.24)  shows  the  disc  compression  in  all 
areas  except  the  region  which  is  not  covered  by  the  bony  endplates.  The 
stress  in  the  z  direction  (fig.  4.25)  also  shows  the  disc  in  compression. 
The  stress  curve  of  the  hoop  stress  shows  to  the  disc  in  compression 

except  in  the  outer  region  of  the  disc.  Once  again,  all  three  stresses 

increase  with  time.  The  z-stress,  shown  in  fig.  4.25,  shows  compression 
throughout  the  disc.  In  addition,  one  can  see  the  stress  again  increas¬ 
ing  with  time.  The  graphs  of  the  hoop  stress,  fig.  4.26  shows  the  disc 

in  compression  in  all  areas  except  the  outermost  area  of  the  disc  where 
it  is  not  covered  by  the  endplate.  Again,  this  is  an  indication  of  the 
shell -like  restraint  of  the  disc's  edge  upon  the  inner  nucleus  as 
discussed  previously. 

4.4.4  Comparison  of  Analytical  Results 

The  comparison  of  the  analytical  results  will  be  done  in  two 
phases:  first  the  healthy  analysis  will  be  compared  to  that  of  the  fully 
herniated  and  second  the  healthy  analysis  will  be  compared  to  that  of  the 
osteo. 

The  comparison  of  the  healthy  analysis  to  the  fully  herniated 
indicates  the  influence  of  the  nucleus  pulposus  on  the  viscoelastic 
material  property.  The  elastic  constant,  E  ,  is  approximately  double 
for  the  healthy,  Eq=  36.5  ^/cm^  compared  to  the  fully  herniated 
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fig.  4.24  Radial  stress  versus  radius  in  the  disc:  of  the 
Os  too  specimen. 
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fig.  4.25  Vertical  stress  versus  radius  in  the  disc  of  the 
Osteo  specimen. 
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Eq=  20.6  ^/cm^.  In  addition,  the  Kelvin  solid  unit  for  the 
herniated  not  only  gave  a  larger  final  displacment  but  the  creep  rate  was 
much  faster  than  the  healthy.  The  stress  curves  also  show  that  the 
increase  in  all  stress  levels  was  greater  for  the  herniated  than  the 
healthy.  From  the  comparison  of  the  material  constants  and  the  stress 
levels,  one  can  conclude  that  the  disc  not  only  acts  as  a  structural 
member  but  also  assist  in  redistributing  the  stress  in  the  specimen  over 
time. 

The  correlation  between  the  osteo  specimen  and  the  healthy 
was  closer  than  that  of  the  herniated  to  the  healthy.  However,  several 
differences  did  exist.  The  osteo  specimen  had  a  significantly  elastic 
spring  constant  than  the  healthy:  Eq  =  36.5  ^D/cm  for  the  healthy  and 
Eq=  56  /cm  for  the  osteo.  Yet  the  Kelvin  solid  unit  in  the 
osteo  crept  at  a  higher  rate  than  the  healthy.  In  comparing  the  stress 
analysis  between  the  two  data  sets,  very  little  difference  is  present  at 
either  the  initial  time  or  for  large  time  values. 
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Chapter  5 


Conclusions 

In  the  final  examination  of  this  problem,  one  can  make  statements 
in  several  areas:  the  experimental  procedure,  viscoelasticity  of  the 
annulus  fibrosis  and  the  nucleus  pulposus,  an  example  of  the  self¬ 
sealing  mechanism,  and  lastly  the  advantages  of  an  axi symmetric 
analysis. 

The  experimental  results  yielded  several  important  conclusions.  By 
cupping  of  the  specimen  before  loading,  it  was  possible,  to  achieve  a 
purely  compressive  displacement  and  avoid  pitch  and  torsion  as 
encountered  by  other  exper imenters.  The  load  was  also  evenly  distri¬ 
buted  across  the  top  edge  of  the  specimen.  By  keeping  the  specimens  in 
a  high  humidity  environment  throughout  the  experimentat ion,  it  was 
possible  to  make  several  load  cycles  on  each  specimen  in  varying  states 
of  herniation  without  inducing  changes  in  material  properties  due  to 
water  loss. 

The  viscoelastic  response  was  present  in  all  specimens.  Yet, 
Kazarian  (11)  has  found  in  certain  cases  that  the  loss  of  the  nucleus 
pulposus  results  in  loss  of  the  viscoelastic  response.  A  possible 
explanation  for  this  would  be  that  a  latent  period  exist  when  the  disc 
continues  to  behave  in  a  viscoelastic  manner  after  loss  of  the  nucleus. 
As  a  result,  one  should  not  discount  the  importance  of  the  disc  in  the 
viscoelastic  response.  The  healthy  disc  also  showed  a  smaller 
displacement  for  the  same  stress  level  than  the  herniated.  A  physical 
interpretation  of  this  would  be  the  healthy  specimen  is  stiffer  and  thus 


undergoes  a  smaller  displacement  for  the  same  load. 

In  the  first  attempt  at  herniation,  the  creep  curve  was  almost 
identical  to  the  healthy.  A  possible  explanation  of  this  was  found  in 
White's  text  (18).  He  stated  that  the  annulus  has  a  self-sealing 
mechanism  to  prevent  loss  of  the  nucleus.  This  study  is  an  excellent 
example  of  the  mechanism.  In  addition,  one  can  now  conclude  that  to 
test  only  the  annulus  one  must  negate  the  nucleus  as  a  structural 
component. 

Finally,  the  advantages  of  an  ax  i symmetric  analysis  give  one  a  feel 
for  the  stress  flow  within  the  specimen.  The  stress  in  the  cortical 
bone  resisting  the  flow  of  the  trabecular  bone  under  compression  shows 
the  cortical  bone  is  behaving  similar  to  a  doubly  curved  anticlastic 
hyperboloid  of  revolution  shell  restraining  the  trabecular  bone. 
Furthermore,  it  appears  that  due  to  the  constraint  imposed  by  the  bony 
endplate  the  outer  region  of  the  disc  is  acting  similar  to  a  doubly 
curved  synclastic  shell  of  revolution  characterized  by  the  tension  in 
the  hoop  direction. 
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Appendix  A 


Experimental  Protocol 

Title:  Creep  Characteristics  of  Normal  and  Herniated  Intervertebral 
Joints  in  the  Lumbar  Region  of  the  Spine. 

Work  Unit: 

Principal  Investigator:  David  R.  Furlong,  2nd  Lt,  AFIT/EN 
Co-Invest i gators:  A.  N.  Palazotto,  Ph.D.,  AFIT/EN 
E.  Paul  France  1  Lt,  BSC 

A.l.  Ob jecti ves:  The  purpose  of  this  investigation  is  to  obtain 
information  on  the  influence  of  the  nucleus  pulposus  on  the  viscoelastic 
properties  of  the  intervertebral  joint. 

A.  2.  Introduct  ion:  Compressive  creep  tests  will  be  performed  on  six 
specimens  from  the  lumbar  section  of  rhesus  monkeys.  The  disc  will  be 
partially  herniated  and  the  creep  test  repeated.  The  disc  will  then  be 
herniated  with  a  trephine  and  the  load  cycle  repeated  for  a  third  time. 
Next,  both  the  experimental  procedure  and  numerical  analysis  will  be 
brought  together  and  further  studied  in  a  Master's  thesis  at  the  Air 
Force  Institute  of  Technology. 

Lumbar  sections  from  six  rhesus  monkeys  at  the  LI-L2  level  will  be 
used  in  the  experiment.  The  general  plan  of  the  protocol  can  be  broken 
down  into  four  parts:  procedure,  equipment  and  instrumentation, 
justification  of  species,  and  results. 

A . 3 . 1 .  Preparation 

1.  The  following  procedure  is  for  the  preparation  of  each 
specimen  and  is  to  be  completed  before  the  morning  of  the  compression 
test. 
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2.  Specimens  will  be  X-rayed  in  both  the  anterior  and  lateral 
position.  These  X-rays  will  be  used  to  establish  colurm  condition  and 
for  anthropomorphic  reference. 

3.  Each  specimen  will  consist  of  the  intervertebral  disc  and 
its  two  adjacent  vertebral  centrums.  All  soft  tissue  and  the  posterior 
processes  will  be  removed  from  the  specimen.  The  total  height  of  the 
test  specimen  will  be  taken  before  and  after  the  posterior  elements  are 
removed.  All  height  measurements,  throughout  the  experimentation, 
should  be  taken  in  both  the  anterior /posterior  and  lateral  directions. 

4.  The  area  of  the  superior  and  inferior  exposed  surfaces  of 
the  specimen  will  be  photographed  in  order  to  obtain  their  areas. 
A. 3. 2.  Phase  I 

1.  The  specimen  will  be  subjected  to  a  creep  test  using  the 
equipment  and  instrumentation  section  of  this  protocol.  A  constant 
compressive  load  of  6.80  Kiloponds  will  be  applied  for  8  hours  followed 
by  16  hours  of  recovery  after  load  removal. 

2.  The  data  during  the  load  cycle  should  be  taken  in  the 
following  manner. 

60  samp les/miri .  0-10  min.  of  load 

6  samples/min.  10-20  min.  of  load 

1  sample/min.  1/3-  8  hrs.  of  load 

60  samples/min.  0-10  min.  of  unload 

6  smaples/miri.  10-20  min.  of  unload 

1  sample/min.  1/3-  8  hrs.  of  unload 

3.  In  addition,  photographs  will  be  taken  throughout  the  load 
portion  of  the  test  and  for  the  first  8  hours  of  the  unload  portion  at 
24  second  intervals. 
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A. 3. 3.  Phase  2 


1.  During  Phase  2,  the  disc  will  be  partially  herniated  and  the 
creep  test  repeated.  At  the  end  of  the  unload  cycle  in  Phase  1,  the 
specimen  will  be  surgically  herniated  with  a  14  gage  hypodermic  syringe 
and  stored  in  a  saline  rich  environment  at  a  cold  temperature.  It  is 
not  to  be  stored  in  saline  solution  because  the  specimen  may  absorb  the 
saline  thus  changing  the  material  properties  of  the  disc.  When  the  disc 
is  herniated,  the  following  measurements  should  be  taken.  Before  and 
after  herniation,  the  total  height  of  the  specimen  is  to  be  taken.  The 
volume  of  fluid  removed  is  also  needed. 

2.  In  the  load  cycle  of  the  partially  herniated  specimen,  data 
will  be  gathered  as  specified  in  Phase  1. 

A. 3.4.  Phase  3 

1.  In  the  last  load  cycle,  the  disc  will  be  fully  herniated  and 
the  creep  test  repeated.  At  the  end  of  the  load  cycle  of  the  previous 
phase,  the  specimen  is  to  be  surgically  herniated  with  a  5mm  trephine 
and  then  stored  in  a  saline  rich  environment.  Before  and  after 
herniation,  the  height  of  the  test  specimen  is  to  be  taken. 

2.  In  the  load  cycle,  the  data  gathered  will  be  in  the  same 
manner  as  the  previous  two  test  cycles. 

A . 3 . 5 .  Phase  4 

1  "fter  the  second  creep  test,  the  specimen  will  be  stored  in 
a  saline  rich  environment  for  at  least  24  hours  before  dissection. 

2.  At  the  time  of  dissection,  the  area  of  the  annulus  fibrosis 
and  bony  endplafe  is  to  be  measured.  In  addition,  the  approximate 
volume  of  the  nucleus  pulposus  will  be  found  through  geometric 


3.  In  the  final  part,  measurements  of  the  vertebral  centrums 
will  be  obtained  as  needed  to  formulate  the  model. 

4.  The  displacement  data  from  each  test  will  be  graphed  using  a 
specially  prepared  plotting  routine. 

5.  The  time  delay  photographs  will  also  be  studied  to  determine 
pitch  and  torsion  of  the  specimen.  Vertebral  unit  will  be  photographed 
to  obtain  cross-sectional  area. 

6.  When  the  specimen  is  cupped,  a  minimum  of  dental  acrylic  is 

to  be  used  to  minimize  shear  in  the  vertebral  bodies.  After  attaching 
the  cups,  the  total  height  of  the  specimen  is  again  to  be  taken. 

7.  Two  pins  will  be  placed  on  the  lateral  surface  of  each 

vertebral  centrum.  The  pins  will  be  used  as  a  reference  fo>~  determining 
pitch  and  torsion  of  the  specimen  during  testing. 

8.  The  specimen  will  then  be  stored  in  a  cold  normal  saline 
rich  environment  for  no  longer  than  two  days  until  tested. 

A. 4.  Equipment  and  Instrumentation:  Intervertebral  joint  creep  tests 
will  be  conducted  using  three  specially  built  chambers.  Each  chamber  is 
designed  to  apply  a  static  vertical  load  to  a  specimen  and  measure  the 
resultant  strain  using  LVDT.  The  load  is  appliedwith  the1  specimen  iri 
an  upright  position  ( super i or / i nf er i or  )  and  can  be  varied  by  the 
addition  or  the  removal  of  lead  weights.  Using  Dana  amplifiers,  the 
LVDT  signal  is  processed  and  becomes  a  calibrated  voltage  output.  The 

output  from  the  three  LVDTs  is  sampled  and  printed  on  a  Silent  700  as  a 

digital  readout.  The  sampling  rate  is  controlled  by  the  operator  who 
has  a  choice  of  between  O.Olsec  and  lO.Omin  rates.  The  LVDT  instru¬ 
mentation  is  calibrated  and  zeroed  before  and  after  each  test  run.  For 
the  24  hr.  duration  of  each  experiment  all  3  test  chambers  are  humid- 


ified  using  a  cool  mist  generator. 

A. 5.  Justification  of  Species:  In  order  to  obtain  meaningful  data 
concerning  the  mechanical  creep  characteristics  of  the  human  spine,  only 
animals  having  a  musculoskeletal  resemblance  to  man  are  appropriate  for 
study.  Because  of  the  great  amount  of  comparative  biomechanical  work 
done  with  the  rhesus  monkey  and  the  baboon  spinal  columns  and  the 
availibility  of  excised  spines  from  cadaver  specimens  of  other 
experiments,  it  is  economically  and  scientifically  prudent  to  utilize 
these  species. 

A. 6.  Results:  The  analysis  and  conclusions  of  this  study  will  be 
presented  as  part  of  a  Master's  thesis  at  the  Air  Force  Institute  of 
Technology. 
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The  experimentation  involved  several  distinct  phases.  First  five  specimens 
were  excised  from  the  L1-L2  level  of  young  adult  rhesus  monkeys  (Macaca  mulatta). 
The  prepared  specimen  consisted  of  the  intervertebral  joint  and  the  two  adjacent 
vertebral  centrums  without  either  the  posterior  elements  or  associated  soft 
tissues.  The  prepared  specimen  was  subjected  to  a  constant  compressive  load  of 
sixty-seven  newtons  for  eight  hours  followed  by  a  relaxation  time  of  sixteen 
hours.  The  specimens  were  then  subjected  to  various  stages  of  herniation  and 
the  creep  tests  were  repeated  after  each  herniation. 

The  analysis  was  also  carried  out  in  several  phases.  First,  the  experimental 
data  was  graphed  using  a  specially  prepared  plotting  routine  for  the  Cyber  750 
machine.  A  mean  displacement  curve  was  then  obtained  for  each  set  of  creep 
tests  by  averaging  the  displacement  vs.  time  curve  of  each  test  phase.  Next, 
all  specimens  were  dissected  and  measured  to  obtain  data  needed  in  formulating 
a  finite  element  model.  In  the  final  phase  of  the  analysis,  a  viscoelastic 
axisymmetric  finite  element  model  was  used  to  quantify  the  experimental  data. 

A  three  parameter  Kelvin  solid  was  employed  in  the  finite  element  solution 
method.  The  result  of  this  paper  will  be  used  to  construct  a  dynamic  model  of 
the  vertebral  column. 
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